Streaming histogram sketching for rapid microbiome analytics

Background: The growth in publically available microbiome data in recent years has yielded an invaluable resource for genomic research, allowing for the design of new studies, augmentation of novel datasets and reanalysis of published works. This vast amount of microbiome data, as well as the widespread proliferation of microbiome research and the looming era of clinical metagenomics, means there is an urgent need to develop analytics that can process huge amounts of data in a short amount of time. To address this need, we propose a new method for the compact representation of microbiome sequencing data using similarity-preserving sketches of streaming k-mer spectra. These sketches allow for dissimilarity estimation, rapid microbiome catalogue searching and classification of microbiome samples in near real time. Results: We apply streaming histogram sketching to microbiome samples as a form of dimensionality reduction, creating a compressed ‘histosketch’ that can efficiently represent microbiome k-mer spectra. Using public microbiome datasets, we show that histosketches can be clustered by sample type using the pairwise Jaccard similarity estimation, consequently allowing for rapid microbiome similarity searches via a locality sensitive hashing indexing scheme. Furthermore, we use a ‘real life’ example to show that histosketches can train machine learning classifiers to accurately label microbiome samples. Specifically, using a collection of 108 novel microbiome samples from a cohort of premature neonates, we trained and tested a random forest classifier that could accurately predict whether the neonate had received antibiotic treatment (97% accuracy, 96% precision) and could subsequently be used to classify microbiome data streams in less than 3 s. Conclusions: Our method offers a new approach to rapidly process microbiome data streams, allowing samples to be rapidly clustered, indexed and classified. We also provide our implementation, Histosketching Using Little K-mers (HULK), which can histosketch a typical 2 GB microbiome in 50 s on a standard laptop using four cores, with the sketch occupying 3000 bytes of disk space

Modulation of host responses by oral commensal bacteria.

Immunomodulatory commensal bacteria are proposed to be essential for maintaining healthy tissues, having multiple roles including priming immune responses to ensure rapid and efficient defences against pathogens. The default state of oral tissues, like the gut, is one of inflammation which may be balanced by regulatory mechanisms and the activities of anti-inflammatory resident bacteria that modulate Toll-like receptor (TLR) signalling or NF-κB activation, or influence the development and activities of immune cells. However, the widespread ability of normal resident organisms to suppress inflammation could impose an unsustainable burden on the immune system and compromise responses to pathogens. Immunosuppressive resident bacteria have been isolated from the mouth and, for example, may constitute 30% of the resident streptococci in plaque or on the tongue. Their roles in oral health and dysbiosis remain to be determined. A wide range of bacterial components and/or products can mediate immunomodulatory activity, raising the possibility of development of alternative strategies for therapy and health promotion using probiotics, prebiotics, or commensal-derived immunomodulatory molecules

The gut microbiome in dogs with congestive heart failure: a pilot study

Compromised gut health and dysbiosis in people with heart failure has received a great deal of attention over the last decade. Whether dogs with heart failure have a similar dysbiosis pattern to what is described in people is currently unknown. We hypothesised that dogs with congestive heart failure have quantifiable dysbiosis compared to healthy dogs that are similar in sex and age. A total of 50 dogs (15 healthy dogs and 35 dogs with congestive heart failure) were prospectively recruited, and their faecal gut microbiome was assessed using 16S rRNA sequencing (Illumina MiSeq platform). There was no significant change in the microbial diversity and richness in dogs with congestive heart failure. However, there was an increase in abundance of Proteobacteria in the congestive heart failure group (p = 0.014), particularly due to an increase in the family Enterobacteriaceae (p = 0.002) and Escherichia coli (p = 0.033). We conclude that dogs with congestive heart failure have dysbiosis, and we show additional trends in our data suggesting that dogs may have a similar pattern to that described in people. The results of this study provide useful preliminary information and raise the possibility that dogs represent a clinically relevant animal model of dysbiosis in people with heart failure

Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences

Profiling phylogenetic marker genes, such as the 16S rRNA gene, is a key tool for studies of microbial communities but does not provide direct evidence of a community’s functional capabilities. Here we describe PICRUSt (Phylogenetic Investigation of Communities by Reconstruction of Unobserved States), a computational approach to predict the functional composition of a metagenome using marker gene data and a database of reference genomes. PICRUSt uses an extended ancestral-state reconstruction algorithm to predict which gene families are present and then combines gene families to estimate the composite metagenome. Using 16S information, PICRUSt recaptures key findings from the Human Microbiome Project and accurately predicts the abundance of gene families in host-associated and environmental communities, with quantifiable uncertainty. Our results demonstrate that phylogeny and function are sufficiently linked that this ‘predictive metagenomic’ approach should provide useful insights into the thousands of uncultivated microbial communities for which only marker gene surveys are currently available

J Acquir Immune Defic Syndr

BackgroundCervical cancer is a major public health problem in resource-limited settings, particularly among HIV-infected women. Given the challenges of cytology-based approaches, the efficiency of new screening programs need to be assessed.SettingCommunity and hospital-based clinics in Gaborone, Botswana.ObjectiveTo determine the feasibility, and efficiency of the \u201cSee and Treat\u201d approach using Visual Inspection Acetic Acid (VIA) and Enhanced Digital Imaging (EDI) for cervical cancer prevention in HIV-infected women.MethodsA two-tier community-based cervical cancer prevention program was implemented. HIV-infected women were screened by nurses at the community using the VIA/EDI approach. Low-grade lesions were treated with cryotherapy on the same visit.ResultsFrom March 2009 through January 2011, 2,175 patients were screened for cervical cancer at our community-based clinic. 253 (11.6%) were found to have low-grade lesions and received same-day cryotherapy. 1,347 (61.9%) women were considered to have a normal examination and 575 (27.3%) were referred for further evaluation and treatment. Of the 1,347 women initially considered to have normal exams, 267 (19.8%) were recalled based on weekly quality control assessments. 210 (78.6%) of the 267 recalled women and 499 (86.8%) of the 575 referred women were seen at the referral clinic. Of these 709 women, 506 (71.4%) required additional treatment. Overall, 264 CIN stage 2 or 3 were identified and treated, and six micro-invasive cancers identified were referred for further management.ConclusionsOur \u201cSee and Treat\u201d cervical cancer prevention program using the VIA/EDI approach is a feasible, high-output and high-efficiency program, worthy of considering as an additional cervical cancer screening method in Botswana, especially for women with limited access to the current cytology-based screening services.20122014-01-08T00:00:00ZP30 AI045008/AI/NIAID NIH HHS/United StatesU2G PS001949/PS/NCHHSTP CDC HHS/United States1U2GPS001949/PHS HHS/United StatesIP30 AI 45008/AI/NIAID NIH HHS/United States22134146PMC388408

Search for rare quark-annihilation decays, B --> Ds(*) Phi

We report on searches for B- --> Ds- Phi and B- --> Ds*- Phi. In the context of the Standard Model, these decays are expected to be highly suppressed since they proceed through annihilation of the b and u-bar quarks in the B- meson. Our results are based on 234 million Upsilon(4S) --> B Bbar decays collected with the BABAR detector at SLAC. We find no evidence for these decays, and we set Bayesian 90% confidence level upper limits on the branching fractions BF(B- --> Ds- Phi) Ds*- Phi)<1.2x10^(-5). These results are consistent with Standard Model expectations.Comment: 8 pages, 3 postscript figues, submitted to Phys. Rev. D (Rapid Communications

Kinetic Modeling of the Assembly, Dynamic Steady State, and Contraction of the FtsZ Ring in Prokaryotic Cytokinesis

Cytokinesis in prokaryotes involves the assembly of a polymeric ring composed of FtsZ protein monomeric units. The Z ring forms at the division plane and is attached to the membrane. After assembly, it maintains a stable yet dynamic steady state. Once induced, the ring contracts and the membrane constricts. In this work, we present a computational deterministic biochemical model exhibiting this behavior. The model is based on biochemical features of FtsZ known from in vitro studies, and it quantitatively reproduces relevant in vitro data. An essential part of the model is a consideration of interfacial reactions involving the cytosol volume, where monomeric FtsZ is dispersed, and the membrane surface in the cell's mid-zone where the ring is assembled. This approach allows the same chemical model to simulate either in vitro or in vivo conditions by adjusting only two geometrical parameters. The model includes minimal reactions, components, and assumptions, yet is able to reproduce sought-after in vivo behavior, including the rapid assembly of the ring via FtsZ-polymerization, the formation of a dynamic steady state in which GTP hydrolysis leads to the exchange of monomeric subunits between cytoplasm and the ring, and finally the induced contraction of the ring. The model gives a quantitative estimate for coupling between the rate of GTP hydrolysis and of FtsZ subunit turnover between the assembled ring and the cytoplasmic pool as observed. Membrane constriction is chemically driven by the strong tendency of GTP-bound FtsZ to self-assembly. The model suggests a possible mechanism of membrane contraction without a motor protein. The portion of the free energy of GTP hydrolysis released in cyclization is indirectly used in this energetically unfavorable process. The model provides a limit to the mechanistic complexity required to mimic ring behavior, and it highlights the importance of parallel in vitro and in vivo modeling

Bacterial community development in experimental gingivitis.

Current knowledge of the microbial composition of dental plaque in early gingivitis is based largely on microscopy and cultural methods, which do not provide a comprehensive description of oral microbial communities. This study used 454-pyrosequencing of the V1-V3 region of 16S rRNA genes (approximately 500 bp), and bacterial culture, to characterize the composition of plaque during the transition from periodontal health to gingivitis. A total of 20 healthy volunteers abstained from oral hygiene for two weeks, allowing plaque to accumulate and gingivitis to develop. Plaque samples were analyzed at baseline, and after one and two weeks. In addition, plaque samples from 20 chronic periodontitis patients were analyzed for cross-sectional comparison to the experimental gingivitis cohort. All of the healthy volunteers developed gingivitis after two weeks. Pyrosequencing yielded a final total of 344,267 sequences after filtering, with a mean length of 354 bases, that were clustered into an average of 299 species-level Operational Taxonomic Units (OTUs) per sample. Principal coordinates analysis (PCoA) plots revealed significant shifts in the bacterial community structure of plaque as gingivitis was induced, and community diversity increased significantly after two weeks. Changes in the relative abundance of OTUs during the transition from health to gingivitis were correlated to bleeding on probing (BoP) scores and resulted in the identification of new health- and gingivitis-associated taxa. Comparison of the healthy volunteers to the periodontitis patients also confirmed the association of a number of putative periodontal pathogens with chronic periodontitis. Taxa associated with gingivitis included Fusobacterium nucleatum subsp. polymorphum, Lachnospiraceae [G-2] sp. HOT100, Lautropia sp. HOTA94, and Prevotella oulorum, whilst Rothia dentocariosa was associated with periodontal health. Further study of these taxa is warranted and may lead to new therapeutic approaches to prevent periodontal disease.BBSRC Industrial Case Studentship ref no. BB/G01714X/1 in collaboration with GlaxoSmithKline

A computational evaluation of over-representation of regulatory motifs in the promoter regions of differentially expressed genes

BACKGROUND: Observed co-expression of a group of genes is frequently attributed to co-regulation by shared transcription factors. This assumption has led to the hypothesis that promoters of co-expressed genes should share common regulatory motifs, which forms the basis for numerous computational tools that search for these motifs. While frequently explored for yeast, the validity of the underlying hypothesis has not been assessed systematically in mammals. This demonstrates the need for a systematic and quantitative evaluation to what degree co-expressed genes share over-represented motifs for mammals. RESULTS: We identified 33 experiments for human and mouse in the ArrayExpress Database where transcription factors were manipulated and which exhibited a significant number of differentially expressed genes. We checked for over-representation of transcription factor binding sites in up- or down-regulated genes using the over-representation analysis tool oPOSSUM. In 25 out of 33 experiments, this procedure identified the binding matrices of the affected transcription factors. We also carried out de novo prediction of regulatory motifs shared by differentially expressed genes. Again, the detected motifs shared significant similarity with the matrices of the affected transcription factors. CONCLUSIONS: Our results support the claim that functional regulatory motifs are over-represented in sets of differentially expressed genes and that they can be detected with computational methods